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The microcanonical ensemble is in important phys- 
ical situations different from the canonical one even in 
the thermodynamic limit. In contrast to the canonical 
ensemble it does not suppress spatially inhomogeneous 
configurations like phase separations. It is shown how 
phase transitions of first order can be defined and clas- 
sified unambiguously for finite systems without the use 
of the thermodynamic limit. It is further shown that 
in the case of the 10-states Potts model as well for 
the liquid-gas transition in Na, K, and Fe the micro- 
canonical transition temperature, latent heat and in- 
terphase surface tension are similar to their bulk values 
for ~ 200 — 3000 particles. For Na and K the number of 
surface atoms keeps approximately constant over most 
of the transition energies because the evaporation of 
monomers is compensated by an increasing number of 
fragments with > 2 atoms (multifragmentation). 

PACS numbers: 05.20.Gg, 05.70Fh, 64.70.Fx, 68.10.Cr 



In the thermodynamic limjv— »«> \^/y =g the mi- 
crocanonical (ME) and the (grand)canonical en- 
sembles (CE) are usually considered to be identi- 
cal. As this is not the case e.g. at phase transi- 
tions of first order and as both ensembles describe 
different physical situations in the case of finite 
systems it is necessary to emphasize here the dif- 
ferences between the two. It is important to re- 
alize that in the microcanonical ensemble phase 
transitions can very well be defined and classified 
for finite systems (sometimes for systems of some 
hundred particles) without the use of the thermo- 
dynamic limit. It is thus possible to define phase 
transitions even in systems which have no thermo- 
dynamic limit at all like systems with unscreened 
forces of long range. The arguments are not en- 
tirely new, however, many discussions with work- 
ers from different physical disciplines showed that 
it is necessary to state these facts clearly. 

The difference between the ME and the CE can 
be seen most easily for the probability P(E/N) 
in the CE at a phase transition of first order 
(T = T tr ) e.g. liquid to gas. Here P(E/N) is a 
bimodal distribution with a peak at the specific en- 
ergy of the liquid Ei/N and another at the specific 
energy of the gas E g /N, As E g /N - Ei/N = q lat 



the specific latent heat, the fluctuations of the total 
energy per particle remain finite even in the ther- 
modynamic limit. Consequently, the canonical is 
different from the microcanonical ensemble. 

The ME is the ensemble which is obtained di- 
rectly from the mechanics of the N-body system. 
The classical partition sum W(E) is the volume of 
the energy shell in the N-body phase space in units 
of (2ttTi) 3N . In the quantum case it is the number 
of N-body states at the energy E. The canonical 
partition sum is obtained from W(E — Ne) by the 
Laplace transform 



Z(T,P,N) = 



W(E, V, N)e- N{e+Pv ^ T dE dV. 

(1) 



As the Laplace transform is a very stable trans- 
formation, the inverse is highly unstable. Un- 
certainties or "reasonable" approximations of the 
canonical Z(T) can imply serious defects in the 
microcanonical W(E). An example was discussed 
for the Bethe nuclear level-density formula in Q. 
In this respect the ME is the fundamental ensem- 
ble. The epigraph on Boltzmann's gravestone is 
the most concise formulation of thermodynamics: 



S = k* InW 



But there is another reason why the ME ap- 
plies for much more physical situations than the 
CE. The microcanonical caloric function T(E/N) 
depends usually only very weakly on the number 
of particles N and reflects already for a couple of 
some 100 particles quite often the bulk properties. 
Moreover, the ME can be applied to inhomoge- 
neous systems like e.g. when there are forces with 
a range longer than the linear dimensions of the 
system. This is the situation in selfgraviting sys- 
tems but also small systems like hot nuclei are tem- 
porarily equilibrated systems under the long range 
Coulomb force. Statistical multifragmentation of 
hot nuclei must be and was treated by microcanon- 
ical statistics from its early days, see The 
CE has the tendency to suppress inhomogeneities 
the larger the number of particles Q , see below. 
It is interesting to notice that Georgii has shown 
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that if one enforces translational invariance of the 
ME, i.e. if one projects out inhomogeneous con- 
figurations (with phase separations) g it becomes 
equivalent to the CE in the thermodynamic limit 
also at phase transitions of first order. 

At phase transitions of first order the system be- 
comes inhomogeneous and fragments into spatially 
separated pieces of "liquid" and "gas" phases. To 
see how the ME handles this situation we discuss 
the Potts model |6|j^] , a generalization of the Ising 
model, defined by the Hamiltonian: 



shaded region whereas c(T) — N(< e 2 > — < 
e > 2 )/T 2 must be positive in the CE. The canon- 
ical ensemble of the bulk jumps over the shaded 
region between the vertical lines at e\ and £3. 
This is the region of coexistence of two phases 
one with ordered spins, the other with disordered 
spins. Here c(e) has two poles and becomes neg- 
ative in between. Notice that the poles are inside 
E\ < e < £3, i.e. the canonical specific heat re- 
mains finite and positive as it should. 
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(2) 



on a two dimensional lattice of N = L * L spins 
here with q = 10 possible components. The sum is 
over pairs of nearest neighbor lattice points only. 
The Potts model can easily be simulated numeri- 
cally and its behavior in the thermodynamic limit 
is known analytically Q . The microcanonical par- 
tition sum over all possible different configurations 
v with the same total energy E is: 



W(E)=J2Se v ,e 



(3) 



Fig.|l|b shows the microcanonical caloric equa- 
tion 1/T = (3(e = E/N), obtained by simulating 
j3 = ds/de in « 12 * 10 5 microcanonical sweeps per 
energy step for a L — 200 lattice. The phase tran- 
sition is manifested by the backbending of f3(e) M. 
The "Maxwell" line defines the transition temper- 
ature and the two shaded areas below(above) are 
just the entropy loss (gain) As sur f when one sep- 
arates(joins) the two phases by an interphase sur- 
face. Spins on the surface have a reduced freedom 
and consequently the entropy of the lattice is re- 
duced per surface spin by As sur f * L = a SU rf jTtr 
the specific surface tension. The length of the 
"Maxwell" line £3 — £1 is the latent heat qi a t per 
particle. This S-bending of the microcanonical 
caloric curve (3(e) is a signal of a phase transition 
of first order [pp^j3]|. 

The total entropy per lattice point is given by 
s(e) — J*g (3(e')de (fig.|l|a). In order to visualize 
the anomaly of the entropy a linear function a + be 
(a = 0.17, b = 1.42) was subtracted. The depth 
of the convex intruder is again the loss of entropy 
As sur f when a surface is created to separate the 
two phases. As we use periodic boundary condi- 
tions one needs two cuts to separate the phases. 



The specific heat c(e) 
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shown in fig.|l|c after smoothing the statistical fluc- 
tuations of (3(e). c(e) becomes negative in the 
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FIG. 1. Microcanonical phase transition of first or- 
der in the Potts model, L = 200, q = 10. 

In the CE configurations with different energies 
contribute cx exp[— N * f(e)/T] to the canonical 
partition sum with f(e) = e — Ts(e) the free en- 
ergy per particle. Due to the convex intruder 
in s(e) by As sur f configurations with coexistent 
phases and energies e ~ £2 are suppressed in the 
CE by the exponential factor exp(—NAs sur f) — 
exp(—o~ sur fL/T). This fact is frequently used in 
canonical model simulations of phase transitions 
to obtain the surface tension e.g. |0,E3. 

In Fig-H the microcanonical surface tension is 
shown as function of 1/L. The surface tension de- 
viates for a lattice of L — 12 by only 10% from 
the bulk values. Other parameters deviate even 
less from the bulk value, e.g. the transition tem- 
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perature by less than 2%. Within the ME a phase 
transition can be identified with all its salient pa- 
rameters, the transition temperature, the specific 
latent heat, the interphase surface tension, and the 
change of entropy per particle. It can unambigu- 
ously be identified and classified for finite systems 
from its microcanonical caloric equation of state 
T(E/N). There is no need to refer to the ther- 
modynamic limit. Evidently, the usual signal for 
a phase transition as a singularity in a canonical 
observable as function of the temperature is arti- 
ficial and applies to the canonical ensemble only, 
c.f e.g. /pl|,p|,p|/. This conclusion has far reach- 
ing consequences: One can easily extend the con- 
cept of a phase transition to systems which do not 
have any thermodynamic limit like e.g. selfgravit- 
ing systems or charged systems under unscreened 
Coulomb repulsion like hot nuclei |2| . 

"Specific Surface Entropy" As*L = AS/L 
of 2dim 10 -states Potts model 
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FIG. 2. Surface tension as function of 1/L. §0,pj 



In order to test the relation of the area under the 
backbend of f3(e) to the empirical bulk surface ten- 
sion for a realistic system, we simulated the phase 
transition for evaporation and multifragmentation 
of some hundreds Na, K, and Fe atoms under 
normal pressure by our microcanonical Metropolis 
Monte Carlo sampling MMMC. The phase transi- 
tion shows up as a clear anomaly (backbend) of the 
microcanonical caloric curve T(e) c.f. fig||. It is 
easy to determine the transition temperature T tr , 
the latent specific heat qi a t and the surface entropy 
per particle As sur t. The transition is clearly iden- 
tified to be of first order. In order to calculate 



the surface tension one has to determine the total 
surface area of all fragments. This is somewhat 
difficult as the size of the fragments drops strongly 
with rising excitation energy because of an increas- 
ing evaporation of monomers. In the fragmenta- 
tion of Na, K clusters the total surface area keeps 
nevertheless roughly constant and close to N ' 
in the first part of the transition in E\ < £ < £2- 
Here the reducing size nii of the fragments is com- 
pensated by an increasing number of fragments 
Nfr — J2Nmi>2- Therefore, the surface tension 
may be estimated by formula (|j), where the aver- 
age is taken over the energy interval £1 < £ < £2- 
As iron decays by pure monomer evaporation be- 
cause of its much larger surface tension, there is 
only one evaporation residue with a steadily di- 
minishing surface area and the use of formula (J4J) 
to determine the surface tension is problematic. 
The values of the so determined surface tension 
approach systematically the bulk values, tabic |. 
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FIG. 3. Caloric curve (thick), number Nf r of frag- 
ments with rrii > 2, and "number" of surface atoms 
^2( m i Ni) of microcanonical systems with 200 and 
1000 sodium atoms at normal pressure. 

Conclusion: Whether an interacting many-body 
system has a phase transition is not a property of 
an infinitely large system in the thermodynamic 
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limit. One can see and classify phase transitions 
in finite systems by the form of the microcanonical 
caloric equation of state T{E/N). This opens the 
possibility to define and discuss phase transitions 
also in systems which are not "thermodynamically 
stable" in the sense of van Hove @. MMMC Jl7]|| 
allows to calculate T(E/N) for some 100 to 1000 
particles. For our examples the transition repre- 
sents quite well the properties of the bulk. Config- 
urations with coexistent, separated phases are well 
represented in theMB whereas in the CE they are 
suppressed by a factor exp(—N 2 / 3 a sur f/T tr ). 





iVo 


200 


1000 


3000 


bulk 


Na 


T tr [K] 


816 


866 




1156 


qiat [eV] 


0.791 


0.871 




0.923 


Sboil 


11.25 


11.67 




9.267 


surf 


0.55 


0.56 








39.94 


98.53 




oo 


V/Ttr 


2.75 


5.68 




6.34 


K 


Ttr [K] 


697 


767 




1033 


qiat [eV] 


0.62 


0.7 




0.80 


Sboil 


10.35 


10.59 




8.99 


surf 


0.65 


0.65 






JV e// 


32.52 


92.01 




oo 


u/Ttr 


3.99 


7.06 




7.40 


Fe 


Ttr [K] 


2600 


2910 


2971 


3158 


qiat [eV] 


2.77 


3.18 


3.34 


3.55 


Sboil 


12.38 


12.68 


13.1 


13.04 


A.S surf 


0.75 


0.58 


0.77 




JV e/.f 


28.34 


61.41 


156.94 


oo 


u/Ttr 


5.29 


9.44 


14.72 


16.57 



TABLE I. Parameters of the liquid - gas transition 
at 1 atm. in a microcanonical system of No interacting 
atoms and in the bulk. The bulk values of a are cal- 
culated by a = 47rr^ s (a a — b s Tt r ) with a s consistent 
to the value used in this calculation for the binding 
energies of the cluster. It was determined by |ll| from 
experimental binding energies of different clusters. It 
is larger than the bulk values of Jisfl . As b s for clus- 
ters is not known and not used in this calculation it 
was taken from JTjfl . For multifragmenting Na and K 
the calculations are considerably more time consuming 
prohibiting the calculation for Ao = 3000. 
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